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Abstract. Classically oscillating massive fields can be used as "standard clocks" in the 
primordial universe. They generate features in primordial density perturbations that directly 
record the scale factor evolution a(t). Detecting and measuring these "fingerprint" signals is 
challenging but would provide a direct evidence for a specific primordial universe paradigm. 
In this paper, such a search is performed for the power spectrum of the Cosmic Microwave 
CN , Background (CMB) anisotropies using the WMAP7 data. Although a good fit to the data 

privileges a scale around k = 0.01 Mpc -1 , we do not find statistical significance for, neither 
against, the presence of any feature. We then forecast the expected constraints a Planck-like 
CMB experiment can impose on the fingerprint parameters by using Markov-Chain-Monte- 
Carlo (MCMC) methods on mock data. We exhibit a high sensitivity zone for wavenumbers 
ranging from 0.01 Mpc -1 to 0.1 Mpc -1 in which fingerprints show up first on the posterior 
■ probability distribution of the wavenumber at which they occur, and then on the modulation 

frequency. Within the sensitivity zone, we show that the inflationary paradigm can be inferred 
from a single feature generating at least a 20% modulation of the primordial power spectrum. 
This minimal value sensitively depends on the modulation frequency. 
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1 Introduction 

Experimentally distinguishing the primordial universe paradigms that lead to the Big Bang 
model is an outstanding challenge in modern astrophysics and cosmology. The leading can- 
didate is inflation [1-3]. While we are still gathering experimental evidences to distinguish 
the inflationary paradigm from other possible alternatives, within the last 15 years we unex- 
pectedly discovered that our late-time universe is actually inflating. One important reason 
this discovery is so convincing is that we are able to directly measure the scale factor of the 
universe a as a function of time t. Using the type la supernovae as "standard candles" [4, 5] 
the measurement of the magnitudes versus the redshifts of the stars directly tells us a(t) - 
the definition of the evolutionary paradigm of the universe. In contrast, the information that 
we have obtained so far from the primordial density perturbations, such as the approximate 
scale-invariant power spectrum, are convoluted consequences of the scale factor evolution, 
and this is a primary reason for possible degeneracies. So, can we directly measure the scale 
factor as a function of time for the primordial universe? 

It has been recently proposed that we may look for "standard clocks" [6, 7]. Such clocks 
should have a known time-dependence and leave their "ticks" in terms of features in the pri- 
mordial density perturbations. They should exist as general as possible in all paradigms and 
leave identifiable characteristics in the density perturbations. Good candidates are classically 
vibrating massive fields. 

By massive, we mean the masses of these fields are much larger than the event-horizon 
energy-scale 1 during the primordial epoch. Such fields are abundant in any primordial uni- 
verse models, for instance in terms of stabilized moduli. The low energy effective field tra- 
jectory, driving the evolution of the universe, is running in the valleys determined by these 
massive fields. All these fields span a multi-field space with very large dimensions. Gener- 
ically, one expects the low energy trajectory to turn from time to time in this multi-field 
space, and, depending on the sharpness of the turns, some massive fields orthogonal to the 

1 For non-inflationary cases, the event horizon energy-scale may not be close to the Hubble scale. 
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adiabatic field trajectory may get excited and oscillate for a while. Such processes have vari- 
ety of manifestations in the low energy theory, appearing as turning, sharp features, particle 
interactions and etc. The induced oscillations typically have small amplitudes. For most 
purposes, they can be safely averaged out or treated as some small side-effects. However, 
these side-effects contain very valuable information. How the massive fields oscillate in a 
time-dependent background can be computed precisely and have several very distinguishable 
features. These oscillations generate "cosmological ticks" that can be used as the above 
mentioned standard clocks. 

The next questions are, how large observational effects can be induced by these small 
oscillations, and how model-independently can we make theoretical predictions? It is shown 
in [6, 7] that three universal properties nicely fit into each other for our purpose. Firstly, 
these oscillations imprint standard clocks in various cosmological parameters in terms of small 
oscillating components; and these parameters appear as couplings in the correlation functions. 
Secondly, these oscillations affect the density perturbations through the universal Bunch- 
Davies vacuum of the quantum fluctuations, instead of their highly model-dependent event- 
horizon scale and super-event-horizon evolution; and this makes general analyses possible 
for different paradigms. Lastly, the sub-event-horizon scale is precisely the place where the 
strong resonance mechanism takes effect; and this greatly enhances the observability of such 
signals for certain parameter space, even if the vibrating field couples to the curvaton 2 only 
through gravity. 

These signals show up as fine-structures in the density perturbations. For the power 
spectrum, they appear as oscillatory corrections to the leading scale-invariant shape. The 
fraction AP^/P^ is typically given by [6, 7] 



Given p, whether the cosmological phase is expanding or contracting is determined by 
the requirement that the quantum fluctuations should exit the event-horizon during this 
epoch. For example, \p\ > 1 corresponds to inflation, in which p > 1 has slowly decreasing 
H (with t > 0) and p < — 1 has slowly increasing H (with t<0);p = 2/3is the matter 
contraction phase; < p <C 1 is the ekpyrotic (slowly contracting) phase; and — 1 <C p < 
describes a slowly expanding phase. For recent reviews on these alternative scenarios, see 
Refs. [8-10]. The parameter m a is the mass of the massive field, k T is the first resonant 
mode excited by the oscillation, and H* is the Hubble parameter at that moment. For the 
expanding background, p > 1 and p < 0, lower A;-modes resonate earlier and the above 
formula applies to 2k > k T ; for the contracting background, < p < 1, larger A;- modes 
resonate earlier and it applies to 2k < k T . The corresponding patterns are illustrated in 



2 In this context, we define the term "curvaton" as the field that sources the leading scale-invariant power 
spectrum. It is one of the field directions in the low energy effective field space mentioned above. We have 
avoided using the term "inflaton" because we are not just considering the inflationary paradigm, also because 
even for inflation the inflaton may not be the dominant source for the density perturbations. 





The parameter we are interested in is p - the index of the fingerprint of the primordial 
universe paradigm - defined through the scale factor evolution as 




(1.2) 



Fig. 1. 
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Figure 1. "Fingerprints" of different paradigms induced by the "standard clocks" in the power 
spectrum AP^/P^. From top to bottom: p = 10 (inflation), p = 2/3 (matter contraction), p = 0.2 
(slow contraction, or Ekpyrosis), p = —0.2 (slow expansion). The Green/light spectra are generated 
by a sharp feature at ko = 100 and have sinusoidal running; these features are not associated with 
"clocks" and their behavior is similar for different paradigms. The blue/dark spectra correspond to 
two massive fields (m CT = 25, 60) excited by this sharp feature and have resonant running; they contain 
the "standard clocks" and can be used to measure a(t) directly. The overall amplitudes of different 
signals have been normalized for clarity. 



The resonant running refers to the scale-dependence in the sine function in Eq. (1.1). 
As we can see, the dependence of this running pattern on the comoving momentum k di- 
rectly measures the time dependence of the scale factor - the two power-law functions are 
inverse to each other. This resonant running behavior is a very robust signature for different 
paradigms because the "zeros" in these oscillations cannot be changed by some potentially 
existing curvaton-isocurvaton couplings during the multi-field evolution. These zeros are 
the cosmological ticks imprinted by the standard clocks. For examples, as we can see from 
Fig. 1, for an expanding background such as inflation, the spacings between the ticks increase 
in a specific way, while for a contracting background such as the matter contraction, they 
decrease. In addition, the paradigms with fast-evolving scale factors, such as the inflation 
and matter contraction, have much more ticks than those with slowly-evolving scale factors, 
such as ekpyrosis. The overall envelop behavior of these oscillatory signals are less robust 
against the multi-field complexities, but their overall scale-dependent trends are very clear. 
Therefore by identifying the running patterns of these signals, determined by the parameter 
p, we measure the fingerprint of a specific primordial universe paradigm. 

Specializing to inflation, there are many types of feature models that can have interesting 
phenomenological consequences, and they give different types of oscillatory signals in density 
perturbations. But it is worth to emphasize that most of them cannot be used as the 
distinguisher between the inflation and the alternative paradigms. 1) The oscillations of 
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massive fields we are interested in are induced by certain sharp features. A sharp feature in 
itself generates a type of oscillatory signals [11-21]. However, these signals are different from 
those generated by the subsequent oscillations of the massive fields. Since a sharp feature has 
only one "click" , it does not provide a "clock" and the running pattern are universally varying 
as 3 ~ sin(2k/ko + phase) for all kinds of background evolution. 2) Some inflation models 
predict small and repeated structures that can generate resonant feature patterns [22-29] . If 
these features are periodic, they generate the same type of oscillatory behavior as from the 
standard clocks, but the scale dependence of the envelop amplitudes are different. However, 
such "clocks" are not standard. Instead, we are probing the properties of these clocks by 
assuming the inflationary background. Phenomenology of non-Bunch-Davies vacuum states 
in inflationary background [30-41] (or non-inflationary background [42, 43]) also belong to 
this class - the non-standard clock is now determined by the property of the new physics 
scale. 3) The de Sitter inflationary phenomenology of oscillating massive fields [44-46] can 
be easily recovered by taking the large p limit in (1.1), in which the power law dependence 
p(2k) 1 ^ p — p becomes the logarithmic dependence ln(2fc). 

In this paper we study how the signals generated by these standard clocks may be 
observed in the cosmic microwave background (CMB) data. While we expect that such 
signals exist generically, their observability depends on the parameter space of the models, the 
state-of-art experimental technologies and data analyses methods. This situation is similar 
to that for the tensor mode, but the signals are of completely different types - the features we 
are looking for are fine-structures in the scalar density perturbations. In addition, the theory 
only predicts the patterns of the fingerprints, but not the index, locations, frequencies, and 
amplitudes. Therefore, a search in all the parameter space is needed. Using CMB data from 
the Wilkinson Microwave Anisotropies Probe (WMAP) satellite, and Planck-like mock data, 
we would like to find out which parts of the parameter spaces are potentially observable and 
can be used to determine the range of the fingerprint index p. 

The paper is organized as follow. Using MCMC methods, we discuss in Sec. 2, the 
constraints set by the WMAP satellite on the eventual presence of signals given by Eq. (1.1). 
Although we do not find evidence for primordial fingerprints in the current data, there is a 
peculiar scale around 10 -2 Mpc" 1 at which such a power spectrum modulation enhances the 
likelihood. 

We then present forecasts for a typical Planck-like CMB experiment. As the likelihood 
for superimposed oscillations in known to be non-Gaussian [38, 39], we use MCMC methods 
on generated mock data for various values of the resonance parameters. This is the subject of 
Sec. 3. We show that there exists a "high sensitivity zone" for 10 -2 Mpc -1 < k < 1CP 1 Mpc -1 
in which primordial fingerprints can be detected 4 . Moreover, in presence of a low signal- 
to-noise feature in the data, the first parameter to be estimated would precisely be the 
wavenumber at which the modulation occurs 5 , whereas frequency, amplitude, phase and p 

3 Here fco parameterizes the location of the sharp feature signal, and k T parameterizes the location of the 
resonant fingerprint signal. The two parameters are related, see Refs. [6, 7]. For p > 0, which include the 
inflation and fast-contraction paradigms, the locations of these two types of signals are hierarchically separated 
(see Fig. 1). 

4 In terms of the multipole moments, using the approximate relation £ ~ kd,A with d,A — 13.7 Gpc [47], the 
sensitivity zone corresponds to 200 < I < 1300. 

5 In [6, 7] and here, the massive field is excited instantly by a sharp feature and decoupled from the 
curvaton afterward (i.e. coupled only gravitationally) . This is the reason there are sharp edges near the 
starting wavenumber fcf in (1.1) and Fig. 1. This makes kf easier to be detected. If the (model-dependent) 
excitation and decoupling process happens more gradually, we expect some smoothing-out effect around the 
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would remain unconstrained. The sensitivity of Planck-like data allows a lower bound on 
the parameter \p\ > 1, the sufficient condition for inflation, to be inferred for amplitudes 
down to max(AP^/P^) ~ 20%. This threshold holds for low frequency signals whereas a 
full reconstruction of the precise value of p requires larger amplitudes, typically greater than 
50%. We also discuss how the reconstruction is affected by the frequency and expansion 
parameter of the underlying signal. We shall conclude in Sec. 4. 



2 WMAP7 data 



In this section, we consider the WMAP seven years data [48-50] for a flat ACDM cosmological 
model with standard parameters, i.e. the density of baryon $7b> of cold dark matter Od m , 
the optical depth r and the Hubble parameter today Hq, or equivalently, 6 (which measures 
the angular size of the sound horizon at last scattering [51]). Concerning the primordial 
power spectrum, we consider only a scalar power-law power spectrum having a "fingerprints" 
modulation as in Eq. (1.1). Defining P c = A; 3 /(2vr 2 ) | C 2 1 , 

we have 



V C = A S 



i + A 



k \ -3+5/(2p) 



sm < (o- 



P 



p-1 



p 



(2.1) 



The parameters ^4 S and n s are the usual amplitude and spectral index while the pivot scale 
has been fixed to its fiducial value fc* = 0.05 Mpc -1 . We have moreover rescaled some of the 
primordial parameters to reduce eventual degeneracies during the data analysis. The scale at 
which the features is observed is fef = fe r /2 and the frequency is now encoded in u = 2m a /H*. 
We have also rescaled the phase compared to Eq. (1.1) as 



tp = + 2 



m CT p" 



(2.2) 



R*P-V 

such that the large p limit does not produce spurious correlations between ip and uj. 
2.1 Parameter space 

In presence of rapid oscillations, the derivation of the temperature and polarization angular 
power spectra require some care. As discussed in Ref. [37], the CMB transfer functions and 
the line of sight integrals have to be estimated with a high precision to avoid under-sampling 
of the oscillations. For this purpose, we have used a modified version of the publicly available 
CAMB code 6 [51]. The price to pay for accuracy is an unacceptable increase of the computation 
time preventing any efficient exploration of the complete parameter space. However, provided 
the modulation amplitude in the Cg remains small enough, and the frequency larger than the 
typical acoustic oscillations, the parameter space associated with the primordial parameter 
remains weakly correlated with the usual cosmological parameters [41, 52]. For this reason, 
and following Ref. [39], we freeze the standard cosmological parameters to their best fit values, 
obtained from a featureless power spectrum, and explore only the space of primordial param- 
eters. The technical advantage being that the CMB transfer functions are computed only 
once. Concerning the WMAP data analysis, we have used the publicly available MCMC code 
CosmoMC [53], together with the likelihood provided by the WMAP team [50] and coupled to 



edges. The universal property of the fingerprints we emphasized here will show up as soon as the massive 
field and the curvaton can be treated as approximately decoupled. 

6 Available at: http : //theory .physics .unige . ch/~ringeval/upload/patches/f eatures 
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Figure 2. Marginalized posterior probability distributions for the primordial parameters (solid), and 
mean likelihood (dotted), associated with the WMAP7 data (at fixed cosmological parameters). There 
is no bound on any fingerprint parameter, i.e. there is no evidence for, neither against, the presence 
of primordial features. Notice however the existence of two favoured scales around k{ ~ 10 _2 Mpc _1 
(but with no statistical significance). 



our modified CAMB code. The sampling is thus performed over the primordial parameters, i.e. 
ln(.A s ), n s , A u , to, p, log(fef) and i/j, starting with a flat prior distribution. The marginalized 
posterior probability distributions are presented in the next section. 

2.2 Constraints on the fingerprint parameters 

As above-mentioned, the standard cosmological parameters have been fixed to their fiducial 
values, i.e. O b /i 2 = 0.02286, n dm h 2 = 0.115, 9 = 1.044 and r = 0.088 (h = 0.71). The MCMC 
exploration has been stopped according to the R — 1 statistics implemented in CosmoMC [53, 
54], i.e. when the estimated variance between different chains does not exceed 1%. This 
number gives the typical error on the resulting posteriors and have been reached for a number 
of samples around 300000. 

In Fig. 2, we have plotted the marginalized posterior probability distributions for the 
primordial parameters (solid curve). Both A s and n s are well constrained, as expected, and 
centered at the best fit values found by considering only a primordial power law spectrum [48] . 
The variance is somehow reduced, but this is the result of having fixed the cosmological pa- 
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rameters. On the other hand, all the other parameters exhibit flat distributions and therefore 
remain unconstrained. The dotted curves in this figure show the mean likelihood [53], which 
typically traces the location of the good fits. As discussed at length in Ref. [39, 55], the 
two distributions do not match when those better fits require some amount of fine-tuning 
between the parameters. Namely, they are located in small and correlated regions of the 
parameter space such that their integrated weight on the marginalized probability remains 
small. In Fig. 2, this situation appears for multiple values of the parameters. For instance, 
the distribution for n s exhibit various "good fits" whereas the overall probability remains of 
Gaussian shape. The weight of each good fit can be assessed by the small deviation induced 
around the Gaussian distribution. Along these lines, the distribution of u clearly exhibits two 
better fits for frequencies around 50 and 300, but the corresponding marginalized distribution 
does not peak by more than one sigma in these locations. In fact, these two frequencies are 
associated with the two favoured scales visible in the distribution of kf, which although not 
statistically significant, are clearly distinguishable. As we will see in the following, the scale 
kf is actually the most sensitive parameter to the actual presence of a feature having a low 
signal-to-noise ratio. 

The overall probability of having detected primordial fingerprints in the WMAP7 data 
is given by the posterior of A u . Indeed, being marginalized over all the other parameters, 
it gives the probability density distribution of having a resonance of amplitude A u , for all 
frequency, phase, location and power p. As one can check in Fig. 2, the A u posterior is 
mostly flat meaning that it is completely unconstrained. We conclude that there is no 
evidence for primordial fingerprints in the WMAP data. On the other hand, high values 
of A u are not disfavoured such that there is also no evidence against (up to the Occam's 
razor favouring the simplest model [56]). This situation is in fact different than unlocalized 
primordial oscillations, such as those coming from non-standard vacuum initial conditions. 
In that case, large amplitudes are disfavoured because the oscillations are spread over all 
multipoles and, if too large, become incompatible with the observed signal [39, 52, 55]. 
Finally, in Fig. 2, one can notice that the posterior of p is slightly reduced around small 
positive values. These models correspond to fast contraction, which, compared to the other 
expansion paradigms, have tendency to produce more widely spread oscillations in the Cg. 
As a result, too large amplitudes are not admissible and this region ends up being slightly 
disfavoured after marginalization. These models will be discussed in more details in Sec. 3.4. 

In view of these results, it is instructive to discuss how much future CMB data can 
constrain the typical signal associated with features. In particular, do the peaks in the kf 
distribution could be interpreted as hints of primordial fingerprints? In the next section, we 
present forecasts for an ideal Planck-like CMB experiment, using similar MCMC methods 
on mock data. We will see that kf is indeed the most sensitive parameter to an underlying 
localized modulation 7 . We will find out the minimum amplitude A u detectable for different 
fingerprints as we vary the frequency u and fingerprint index p. 

3 Planck-like CMB data 

Given a fiducial set of cosmological parameters {A a }, the easier method to forecast their 
expected errors would be to use the Fisher matrix formalism which merely consists as ap- 
proximating the likelihood by a Gaussian around the best fit location [57]. However, as 

7 Let us notice however that such a modulation, in a realistic experiment, may also be generated by some 
residual colored noise. 
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1 1 

Figure 3. Temperature angular power spectrum, plus noise, for an inflationary feature (left panel) 
located at fcf = 0.05 Mpc -1 , A M = 0.5, Cj = 600, p = 8 (inflation) and '0 = 0. The right panel shows 
a feature at the same location but during a fast contracting era having fcf = 0.05 Mpc -1 , A u = 0.4, 
Cj = 200, p = 2/3 (matter contraction) and ■0 = (see Sec. 3.4). 



discussed in Refs. [38, 39, 55] the likelihood associated with superimposed oscillations in the 
Ci is non-Gaussian and can be multi-valued such that Fisher matrix method can only be ap- 
plied for a high signal-to-noise detection. For this reason, we prefer in the following a MCMC 
exploration based on mock data containing various fingerprint oscillatory patterns [38]. 

3.1 Mock power spectrum and likelihood 

Forecasts can be made through an MCMC exploration of the parameter space provided one 
specifies the likelihood. The mock data {Ci} are assumed to be associated with a set of 
fiducial cosmological and primordial parameters {A a } and one needs to specify the likelihood 
of the theoretical {Ci(X a )} tested. For a full sky analysis, assuming isotropic white noise for 
each pixel and Gaussian statistics, one can show that the sampling distribution followed by 
the Ci is a Wishart distribution [58-61]. Including polarization yields, up to a normalization 
constant [62] 

-2\nC[p t \D^\ =5^(2^ + 1) Tr (t) t ■ DJ 1 ^ - In D e ■ DJ 1 -3 . (3.1) 

The matrix T)i = {Df Y } where X, Y stand for temperature and polarization variables, T, 
E and B. The spectra Df Y include a white noise component compared to the angular power 
spectra Cf Y and are defined by 

nXY /iXY , ^noise fo o\ 

1 = £ {Bf) 2 ' ^ ' 

where Bf is the beam response. For a Planck-like experiment, we have chosen a Gaussian 
beam with a full width at half maximum (fwhm) of 7'. The noise power for temperature 
has been set to Cj oisc = 2 x 10~ 4 /iK 2 and for polarization to C^ oisc ~ C^ oise ~ 2C^ ise [63]. 
From our modified version of the CAMB code, and the above noise specification, the Df Y can 
be computed for any input value of the cosmological and primordial parameters {A a }. Two 
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Figure 4. Marginalized posterior probability distributions for the primordial parameters (solid), and 
mean likelihood (dotted), associated with the mock Planck-like data when there is not any feature 
present. Notice the posterior of fcf which inversely traces the feature sensitivity region. 



examples of the temperature angular power spectrum are represented in Fig. 3. Note that 
the oscillation amplitudes are greatly reduced in the temperature angular power spectrum 
comparing to those in the primordial power spectrum. 

Using MCMC sampling with the likelihood of Eq. (3.1) allows to extract the posterior 
probability distribution for each "measured" parameter A a given the fiducial model {A a }. 
For the same reasons discussed in Sec. 2, we have fixed the cosmological parameters to their 
best fit value, obtained from a standard power-law primordial power spectrum, and ran 
the MCMC exploration only in the primordial parameter space. The standard cosmological 
fiducial parameters are the same as in Sec. 2, plus n s = 0.97 and ln(10 1( M s ) = 3.166. In the 
following, we consider various fiducial values for the fingerprint parameters and discuss how 
well they can be reconstructed. 

3.2 No feature: the sensitivity zone 

The first fiducial model considered has no feature, i.e. A u = 0. Practically, the MCMC 
chains are run as specified in Sec. 2, with the same convergence criteria, i.e. the chains are 
stopped when the expected error on the marginalized distributions does not exceed a few 
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Figure 5. One- and two-sigma confidence intervals in the planes (log k{ , A u ) and (log kf,u>) and mean 
likelihood (shading) for vanishing fingerprints in the mock data. These domains trace the sensitivity 
region of Planck-like CMB data. For k{ > 10" 1 Mpc -1 or fcf < 10 -2 Mpc -1 , localized resonances are 
hardly detectable. 



percents. In Fig. 4, we have represented the marginalized posterior distributions obtained 
from the MCMC exploration of the primordial parameter space. 

The distribution for the standard primordial parameters, n s and ^4 S , are reconstructed, 
as expected for well constrained parameters. Almost all fingerprint parameters exhibit a 
flat distribution and are unconstrained. The slight deviations between mean likelihood and 
marginalized posteriors for A^, uj and fef are the signatures of correlations. In Fig. 5, we 
have plotted the two-dimensional one- and two-sigma confidence intervals as well as the two- 
dimensional mean likelihood (shading). This plot shows that there is a strongly disfavoured 
region for k{ between 10 -2 Mpc -1 and 10 —1 Mpc~ . In this domain, Planck-like data are 
sensitive to the presence of features and as our fiducial model has no feature, this region is 
disfavoured. Conversely, the other domains are poorly constrained. On smaller scales the 
noise starts to dominate whereas on larger scales any resonance pattern is smoothed out by 
the CMB transfer functions [39, 55]. The one-dimensional distributions of Fig. 4 end up 
being flat because, for any A u (or uj) value, there exists a scale kf for which the model can 
be made undetectable within a Planck-like accuracy CMB experiment. 

3.3 Fingerprints of inflation 

As a starting point, we consider an inflation model with p = 8, having a strong feature 8 
of amplitude A u = 0.5, at a high frequency uj = 600, and located in the middle of the 
detectable zone, i.e. fcf = 0.05 Mpc -1 (tp = 0). Such a model has an angular temperature 
power spectrum represented in Fig. 3 (the polarization spectra are not represented but exhibit 
a similar behaviour). 

Running a MCMC exploration on the primordial parameters against such a model yields 
the marginalized distributions of Fig. 6. The power spectrum amplitude ^4 S is always well- 
constrained and its posterior remains identical to the one of Fig. 4. It is not represented in 
the following. The phase ip is poorly recovered whereas all the other fingerprint parameters 

8 Despite of the large amplitude in the primordial power spectrum, in slow-roll inflation this example 
corresponds to a transfer of ~ 3% of the inflaton kinetic energy to the massive field [6] . 
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Figure 6. Marginalized posterior probability distributions for the primordial parameters (solid), and 
mean likelihood (dotted), for typical Planck-like mock CMB data having a feature with = 0.5, 
u) = 600, i/> = 0, p = 8 and fcf = 0.05 Mpc^ 1 (same as in Fig. 3). Being in the high sensitivity zone, 
all fingerprint parameters are well reconstructed, especially the scale fcf = fc r /2. For such a strong 
signal, the fingerprint parameter p can even be precisely measured, in addition to the conclusion that 
\p\ > 1. 



are well determined. Their standard deviation are AA^ = 0.08, Au = 43, A*S> = 1.2 
and Ap = 2.9 showing that such a feature would indeed allow to probe the expansion rate 
through p. The most sensitive parameter remains however the wavenumber scale as we find 
its standard deviation to be Akf/kf = 4 x 10~ 3 . 

In order to test the sensitivity of Planck-like CMB data with respect to the amplitude, 
we have redone the same analysis for various fiducial amplitudes A u , all the other fiducial 
primordial parameters being unchanged. The marginalized posteriors are represented in 
Fig. 7 and shows that for A u < 0.4, any detection becomes unlikely and impossible for 
A w < 0.3. Let us notice that already for A u = 0.4, amplitude, phase and the expansion 
parameter p are poorly, if not recovered. This suggests that, for a high frequency fingerprint, 
the presence of a strong feature, as the one discussed previously, is crucial for probing the 
p parameter. There exists a very sharp line in terms of the value of A u between the signals 
that can be reconstructed and those cannot. On the other hand, the parameters such as the 
frequency u) and the feature scale fcf still let some imprints, down to A u = 0.3. 

The effect from changing the fiducial frequency in the mock data is typical of any other 
resonant pattern: the CMB transfer function smoothing out rapid oscillations, at constant 
A^, higher frequencies produce a lower signal in the CMB [37]. In Fig. 8, we have plotted 
the posteriors obtained for a fiducial model having Cj = 300 and for various values of the 
amplitude A u . The behaviour is exactly the same as for uj — 600, but the sensitivity to A^ 
is increased because the resonances are less smoothed out by the CMB transfer function. So 
signals with lower amplitudes become slightly more accessible. 
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Figure 7. Marginalized posteriors obtained by reducing the feature amplitude to = 0.4 (red 
dashed) and =0.3 (black solid). The other parameters are the same as in Fig. 6. Already for 
A^ = 0.4, the posterior for A\, gets noisy showing that such a feature is no longer detectable. However, 
frequency and scale are still felt and are the most sensitive parameters. Amplitude, phase and the 
expansion parameter p requires a stronger signal to be inferred (see Fig. 6). 



We have further tested examples with relatively low frequencies. In Fig. 9 we have shown 
the posteriors for an underlying fingerprint with w = 100 and = 0.2. We can see that, for 
low frequency, a fingerprint would be detected by Planck with much lower amplitude. For 
example, the posterior for is sharply peaked at the fiducial value with a standard deviation 
of AA^ = 0.03. While most of the parameters can be reconstructed as before, there is an 
interesting exception of the posterior of p. This parameter can now only be constrained from 
below: p > 5 at 95% of confidence. The loss of accuracy on p also affects the determination of 
cj as both parameters are degenerated (see the two-dimensional posterior in Fig. 9). This is 
not difficult to explain. In the large \p\ limit, the fingerprint profile approaches to a unique de 
Sitter limit, so fingerprints with large p values tend to be degenerate. With lower frequency, 
this degeneracy becomes more effective since the "standard clock" is running slower and there 
are fewer "ticks" available to reconstruct the exact parameter p. However for our purpose, 
the exact value of p is not the most important one. As long as we can demonstrate \p\ 3> 1, we 
would be able to identify the inflation as the underlying paradigm. Interestingly, Fig. 9 shows 
that the posterior for p indeed unambiguously indicates that this is an inflationary paradigm, 
although recovering p = 8 is no longer possible. We have also tested a lower frequency signal 
having Cj = 50 (and A w = 0.2, figures not represented). The posteriors are very similar 
to those of Fig. 9, the amplitude and frequency are peaked at their fiducial value. The 
exact p value cannot be reconstructed as well, and now the degeneracies between (uj,p) are 
extended to negative p values (including the other inflationary branch p < —1). The bound 
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Figure 8. Same as Fig. 7 but for a reduced frequency Cj = 300. The sensitivity is slightly increased as 
signals down to A u ~ 0.3 are now detectable in the posterior for A u . Only the constrained fingerprint 
parameters are represented for three fiducial amplitudes A u = 0.3 (dashed blue), A u = 0.2 (solid red) 
and Au =0.1 (thick black). 



for inflation, \p\ > 1, can still be established but by not more than two-sigma as the 99% 
confidence region includes other paradigms. Further lowering the frequency however will not 
help to reduce the minimum detectable amplitude. This is because, with fewer oscillations, 
the fingerprints start to be confused with the acoustic oscillations. They sometimes do not 
even appear as oscillations, but only as deformations of various acoustic peaks. 

In order to explore different possibilities, we have also tested a rather unusual inflation 
model with p = —2 and lu = 300, all the other fiducial parameters being as before. This is 
still inflation in the sense that \p\ > 1, but unusual because the expansion rate is much slower 
than the de Sitter space and the Hubble parameter is increasing. As Eq. (1.1) shows, the 
fc-dependency is such that the "instantaneous" frequency decreases faster than a logarithm 
for k/kf S> 1. From the CMB point of view, it means that the signal is less damped by the 
transfer functions and the oscillatory pattern is spread over larger multipoles than for the 
inflationary paradigm (p ^ 1). In Fig. 10, we have represented the marginalized posteriors 
obtained by a MCMC analysis. Compared to the case p = 8, fingerprint parameters are more 
constrained due to the larger CMB signal, this is particularly clear for the phase tp. Varying 
amplitude and frequency reproduces the same qualitative behaviour discussed before, namely 
one would find that all features disappear for < 0.1. Again, kf remains the most sensitive 
parameter as we find his posterior well peaked down to A w = 0.2. 
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Figure 9. Posteriors obtained with a small inflationary feature at low frequency: ~ 0.2 and 
Cj = 100. Although low frequencies render smaller oscillation amplitude detectable, the p value 
cannot be reconstructed but still indicates an inflationary paradigm. At lower frequencies, for instance 
A u = 50, the posterior of p extends to negative values such that distinguishing the inflationary 
paradigm from another expansion era becomes more difficult (see text). The lower right panel shows 
the one- and two-sigma confidence intervals of two-dimensional posterior in the plane (ui,p): the two 
parameters are correlated. 



3.4 Beyond the inflationary paradigm 

In this section, we discuss the detectability of the fingerprints within alternatives to inflation. 
We consider three categories: slow-contraction models (0 < p < 1), slow-expansion models 
(— 1 <C p < 0) and fast-contraction models (p < 1). Fast-expanding but non-inflationary 
cases (— 1 < p < and p ~ — 1) are similar to the p = — 2 case considered in the previous 
section. 

To test the slowly contracting models, we consider the case p = 0.1. The effective 
resonance frequency being given by cop 2 , superimposed oscillatory patterns end up being of 
observable frequency for u S> 1/p 2 - For this reason, we have considered a fiducial model 
having to = 6000, its temperature angular power spectrum is represented in Fig. 11 for 
= 0.5. The small value of p makes the oscillation pattern very localized around kf 
and this model ends up being unobservable. Even for an unrealistic maximal modulation 
amplitude (A u = 1), the MCMC analysis does not allow the reconstruction of any fingerprint 
parameters. 

The case of slowly expanding models is very close, their fingerprints exhibit similarity to 
the slowly contracting models which would affect their detectability (see the last two figures 
in Fig. 1). 

Faster contracting models, having larger values of p, should not suffer from this problem, 
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Figure 10. Primordial posteriors for an expanding model having p = —2, with = 0.5 and L) = 300. 
The expansion parameter p is recovered at a standard deviation of Ap = 0.15. 




Figure 11. Two types of undetectable features for a Planck-like experiment. Left panel: zoom on 
the second and third acoustic peaks for a slow contracting model p = 0.1, A u = 0.5 and ui = 6000. 
The resonances are before the feature scale fcf = 0.05 Mpc -1 but remain very localized due to the 
small value of p. Right panel: an inflation model, p = 8, having a strong feature A u = 0.8 washed 
out under cosmic variance by the large scale CMB transfer functions. It has fcf = 0.01 Mpc -1 and 
lu = 300 (see Sec. 3.5). 



as for instance a matter contraction with p = 2/3. On the contrary, they exhibit widespread 
oscillations in the primordial power spectra for k < kf but, as seen in Eq. (1.1), their "instan- 
taneous" frequency increases dramatically. Depending on the values of if the frequency 
u> becomes too large, the CMB transfer functions acting as a low-pass filter strongly damp 
the oscillatory tail at large multipoles. As a result, the fingerprints may no longer be visible 
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Figure 12. Expected marginalized posteriors for a feature generated during a fast contracting era 
and having A u = 0.4, Cj = 200, p = 2/3 (the temperature power spectrum is plotted in Fig. 3). 



above some £ value, which could be lower than if, thereby rendering parameter reconstruc- 
tion difficult. These models can therefore be "visible" only for not too large frequencies, but 
this also implies the existence of a few oscillations in the Cg having a stronger amplitude 
than for inflation. This is illustrated in the right panel of Fig. 3 where we have plotted the 
angular power spectrum associated with the feature Co = 200 and A w = 0.4. The resonance 
patterns show up at the top of the first peak whereas the feature scale is located at much 
smaller scales if ~ 700. In fact, such a characteristic explains why those models are a bit 
disfavoured by the WMAP7 data: compared to other early universe paradigms, at same am- 
plitude and frequency, those p values are associated with slightly too large oscillations in the 
Cg. Concerning the Planck forecasts, we have represented in Fig. 12 the posteriors obtained 
from the MCMC analysis. For the same fiducial amplitude A u = 0.4, all parameters are well 
reconstructed and in particular the expansion index p. The above-mentioned sensitivity to 
the model parameters renders smaller amplitudes hardly detectable for Planck-like data. We 
have tested a smaller fiducial amplitude of A u = 0.3 (figures not represented) for which the 
posterior of A u appears very noisy and not clearly peaked. However, and contrary to the 
inflation paradigms, the posterior of p still indicates a contracting model because, even noisy, 
the oscillatory patterns remain typical of a fast contracting background. In that situation, p 
ends up being a more sensitive parameter than the amplitude. For A u < 0.2, all hints for a 
signal are lost and the posteriors are identical to those of Fig. 4. 

3.5 Outside the sensitivity zone 

As one may expect, if the primordial feature occurs at a scale kf which is well outside 
the sensitivity region, it cannot produce a large enough signal in the CMB to be clearly 
distinguished. To check how the transition occurs, we have performed a MCMC exploration 
for the same fiducial model as in Sec. 3.3 (p = 8, Cj = 300) but with a scale kf = 0.01 Mpc -1 . 
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Figure 13. Marginalized posteriors for a feature occurring at the small scale boundary of the sensi- 
tivity zone, kf = 0.1 Mpc~ . Most of the fingerprint parameters are very well reconstructed but not 
p and uj which are strongly degenerated. Their respective one- and two-sigma confidence intervals 
are plotted in the rightmost lower frame. This is the result of the noise amputating some part of the 
oscillatory tail. Notice that the inflationary paradigm can still be inferred as \p\ > 1. 



Even for an amplitude A u = 0.8, the marginalized posteriors for all primordial parameters 
are identical to those of Fig. 4, i.e. as if no feature were present. As explained before, this is 
the result of the strong smoothing by the CMB transfer functions, which are very efficient on 
the largest scales. Such a smoothing is visible on the fiducial power spectra as the resonance 
patterns remain under the cosmic variance (see right panel of Fig. 11). Let us mention that 
not increasing the numerical accuracy for the transfer functions, i.e. using CAMB at its default 
numerical precision, may produce spurious enhanced signals in the CMB [37]. 

At smaller scales, the noise dominates. In Fig. 13, we have represented the posteriors 
obtained when the feature lies at the small scales boundary of the sensitivity domain, i.e. for 
h = 0.1 Mpc -1 (£{ ~ 1400) and for p = 8, A w = 0.5, u = 300. There is a net signal detection, 
the posterior of A w is sharply peaked around the expected fiducial value (A^4 W = 0.04), as well 
for the scale kf (Afef/fcf = 9 x 10 ). However, the frequency probability distribution appears 
to be multi-valued as well as the one for the expansion parameter p. In the same figure, we 
have plotted the one- and two-sigma confidence intervals associated with the two-dimensional 
posterior probability distribution in the plane (u>,p). Both parameters end up being strongly 
correlated as all shaded value provide a good fit to the CMB resonance pattern. This is 
reminiscent with the degeneracy mentioned previously when the frequency of an inflationary 
fingerprint is low. Here, compared to the same feature in the sensitivity zone, part of the 
oscillatory tail is actually truncated as becoming of smaller amplitude than the noise (see 
Fig. 3). As a result, there are less measurable "ticks" in the CMB and all combination of p 
and u producing similar oscillations around kf cannot be distinguished. It is therefore not 
surprising that the posteriors obtained here are similar to the ones of Fig. 9. Let us notice 
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that the data still "see" that this is an inflationary era - the oscillations remain on the right 
side of kf, independently of the truncation, and the posterior for p is non- vanishing only when 
\p\ > 1. Further out of the sensitivity zone, the fingerprint signals remain undetectable. 

4 Conclusions and discussions 

We have studied the detectability of the fingerprints induced by standard clocks in various 
primordial universe paradigms. Concentrating on the power spectrum, we have found that 
they are detectable in Planck-like data provided they lie in the sensitivity zone 10 -2 Mpc -1 < 
k{ < lO -1 Mpc" 1 . The scale kf, and frequency oj are the most sensible parameter to the actual 
presence of a fingerprint in the data in most cases. For inflation, although their posteriors 
exhibit characteristic signatures for an amplitude as small as = 0.1 (10% modulation), 
we have shown that a proper reconstruction of the parameter p requires a higher signal 
A w ^ 0.5, especially at high frequencies (uj > 300). The situation is however improved 
at lower frequency (to ~ 100) as amplitudes > 0.2 still allow to infer the inflationary 
paradigm \p\ > 1. Alternatives to inflation with small \p\ <C 1 have been found to deviate 
from this rule as very small value of p end up being more difficult to detect. As discussed 
in the previous sections, this is essentially due to how widespread the observable oscillatory 
patterns are. In any case, the high sensitivity zone in kf is the region to explore with the 
soon to be released Planck data. 

We have also seen that, for inflation, establishing the bound \p\ > 1 is easier than a full 
reconstruction of the index p. So the main strength of the standard clocks is to break the 
leading degeneracy between the paradigms with \p\ > 1 and \p\ < 1, especially if they predict 
degenerated power spectrum and non-Gaussianities in the absence of the fingerprints. Once 
this degeneracy is broken, the more precise value of p could be inferred from other more 
standard observables such as the spectral index. 

Most likely, hints for a feature could show up while not allowing a full reconstruction with 
the CMB power spectrum only. However, these standard clocks also imprint correlated signals 
in non-Gaussianities [6, 7]. These signals appear as leading order large non-Gaussianities 
instead of small corrections, due to the resonant mechanism [22]. The search for such scale- 
dependent and non-separable non-Gaussianities is a much more difficult task however (see 
e.g. Ref. [64-66] for possible methods). An efficient approach would be to search them first 
in the power spectrum, as we discuss in this paper. If any candidate signals are found, 
the locations of the corresponding non-Gaussianities will be determined, and the subsequent 
search for non-Gaussianities would be considerably narrowed down and provide non-trivial 
checks. 

Finally, let us mention that the matter power spectrum, in the linear regime, is far less 
smoothed than the CMB concerning the transfer of features. However, those resonances are 
completely washed out by the galaxy survey window functions [37]. A possible future work, 
however, may be to discuss how much structure formation may be affected by features in the 
non- linear regime [67]. 
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